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ABSTRACT 

Context. The protoMIRAX hard X-ray imaging telescope is a balloon-borne experiment developed as a pathfinder for the MIRAX 
satellite mission. The experiment consists essentially in a coded-aperture hard X-ray (30-200 keV) imager with a square array (13x13) 
of 2mm-thick planar CZT detectors with a total area of 169 cm^. The total, fully-coded field-of-view is21°x21° and the angular 
resolution is 1°43'. In this paper we describe the protoMIRAX instrument and all the subsystems of its balloon gondola, and we show 
simulated results of the instrument performance. 

Aims. The main objective of protoMIRAX is to carry out imaging spectroscopy of selected bright sources to demonstrate the perfor¬ 
mance of a prototype of the MIRAX hard X-ray imager. 

Methods. Detailed background and imaging simulations have been performed for protoMIRAX balloon flights. The 3-cr sensitivity 
for the 30-200 keV range is ~1.9 x 10"^ photons cm"^ s"^ for an integration time of 8 hs at an atmospheric depth of 2.7 g cm“^ and an 
average zenith angle of 30°. We have developed an attitude control system for the balloon gondola and new data handling and ground 
systems that also include prototypes for the MIRAX satellite. 

Results. We present the results of Monte Carlo simulations of the camera response at balloon altitudes, showing the expected back¬ 
ground level and the detailed sensitivity of protoMIRAX. We also present the results of imaging simulations of the Crab region. 
Conclusions. The results show that protoMIRAX is capable of making spectral and imaging observations of bright hard X-ray source 
flelds. Furthermore, the balloon observations will carry out very important tests and demonstrations of MIRAX hardware and software 
in a near space environment. 
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1. Introduction 


Coded mask experiments have been widely used in astronom¬ 
ical hard X-ray (E > lOkeV) observations of the sky due to 
several factors. First, the technology to develop focussing tele¬ 
scopes above these energies has only recently been mastered and 
implemente d, especially in the NuSTAR observatory ( [Harrison 


jet al.||2Q13] ). The ASTRO-H mission, expected to be launched 
in the near future, will also make use of focussing optics in 
the hard X-ray range t hrough its Hard X-ray Telescope (HXT) 
( jTakahashi et al.|2014| ). While providing much higher sensitivi¬ 
ties and better angular resolution in general, these instruments 
require very long focal lengths and sophisticated control and 
alignment systems. In addition, the highest energy achievable is 
still significantly below 100 keV, even though there are promis¬ 
ing higher-energy (< 6 00 keV) focussing t echniques like Faues 
lenses being developed ( jVirgilli et al.|2013| ) . Coded mask instru¬ 
ments, on the other hand, allow imaging over very wide fields- 
of-view (which can be a significant fraction of the sky) and up 
to energies of several hundred keV. Another factor is that coded 
mask instruments can be compact and easy to implement, mak¬ 
ing them good options for wide-field hard X-ray and low energy 
y-ray monitors of the highly variable and transient source popu¬ 
lations in this energy range. Important examples of coded mask 
satellite instruments that have made signi ficant contri butions to 
astronomy are: SIGMA ( [Paul et al.|I991| ) and ART-P ( jSunyaevj 
jet al.|I990|) instrum ents onboard GRANAT, t he ASM on RXTE 
( Fevine et al.|1996| ), the WFCs on BeppoSAX ( jJager et al.|I997| ), 
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the WFM on HETE-2 ([Ricker et al.||2003j), instruments on the 
INTEGRAF satellite ( [Winkler[1995] ), and BAT/SWIFT ( [Gehrels[ 
[& Swift Team|2004| ). 

In this paper we describe the protoMIRAX balloon experi¬ 
ment, a wide field hard X-ray coded mask imager under devel¬ 
opment at the National Institute for Space Research (INPE) as 
a pathfinder for t he MIRAX (Monito r e Imageador de Raios X) 
satellite mission ( [Braga et al.[[2004| see also AIP Conference 
Proceedings 840, 2006, for a series of papers and information 
about MIRAX). The main instrument of the balloon payload is 
a hard X-ray camera that employs an array of 13x13 CdZnTe 
(“CZT”) detectors as the position-sensitive detector plane. The 
experiment will be launched on a stratospheric balloon over 
Brazil at a latitude of ~ 23° S to carry out imaging spectroscopy 
of selected bright sources. In addition to testing the camera per- 
fomance in a near space environment, the experiment is a testbed 
for new space technologies being developed at INPE. A new atti¬ 
tude control and pointing system for scientific balloon gondolas, 
including two new star trackers, will be tested, as well as novel 
data handling and ground support systems. 


In sectionwe present the scientific objectives of the exper¬ 
iment. In section we describe the X-ray camera and the other 
subsystems that protoMIRAX comprises. In section]^ we present 
the results of Monte Carlo simulations of the camera response at 
balloon altitudes, showing the expected background level and 
the detailed sensitivity of protoMIRAX. A simulated image of 
the Crab region is shown in section We present our conclu¬ 
sions in section |6l 
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2. Scientific Objectives 

The protoMIRAX experiment will perform imaging spec¬ 
troscopy with good energy resolution of bright Galactic hard 
X-ray sources such as the Crab Nebula (and its pulsar) and the 
Galactic Center complex. 

In X-ray astronomy, the Crab has been considered as a stan- 
dard candle due to its large and nearly constant flux at Earth (see 
IBiihler & Blandford| ( |2014| ) and |Madsen et al.| ( r20I5| ) for recent 
reviews). At X-ray energies above ~30keV, the Crab is generally 
the strongest persistent source in the sky and has a diameter of 
~1 arcmin. The Crab has experienced some recent flares in the 
GeV range ( [Tavani et al. |2Q 11 1 ) which were simultaneous or just 


previous to a parti cularly active period in hard X-rays ( Wilson- 
Hodge et al.|20II| ). Those surprising variations are of the order a 


few percent on the 50-100 keV energy range. ProtoMIRAX will 
be able to measure the Crab spectrum from 30 to 200 keV. These 
observations will also be used for flux calibrations and imaging 
demonstrations. 

The Galactic Center (GC) region is very favorable for pro¬ 
toMIRAX observations due to its declination of d-29°, since 

the balloon will fly over a region at latitudes of ~ -23°. One 
of the most interesting sources is IE 1740.7-2942, the brightest 
and hardest persistent X-ray source within a few degrees of the 
GC. Due to a similar h ard X-ray spectrum and comparable lumi¬ 
nosity to Cygnus X-1 ( [Liang & Nolan||1984| ), IE 1740.7-2942 
is classifled as a black hole candidate. Because of the two-sided 
radio jets associate d to the source, the object was dubbed the 
first “microquasar” ( Mirabel et al.|1992) , an X-ray binary whose 
behavior mimics quasars on a much smaller scale. Since the 
GC direction has extremely high extinction, a counterpart has 
not yet been identified despite deep searches made in both the 
optical and IR bands. The source was recently studied by our 
group from soft to hard X-rays up to 200 keV ( [Castro et al 


[2014[ ), showing spectra that can be very well modeled by ther¬ 
mal Comptonization of soft X-ray photons. We will image the 
Galactic Centre region with protoMIRAX with IE 1740.7-2942 
in the centre of the fleld-of-view. Our observations in the hard 
X-ray band could be useful to measure flux and spectral param¬ 
eters, especially because the source is most of the time in the 
hard state. 

Another interesting source that will be in the same fleld- 
of-view is GRS 1758-258, also a microquasar, and one of 
the brightest X-ray sources near the GC at energies greater 
than 50 keV ( Sunyaev et al.[p^91a[ ). Like its sister source 
IE 1740.7-2942, GRS 1758-258 spends most of the time in the 
hard state. Its hard X-ray spectra and varia bility are also similar 
to that of CygX-1 ( Kuznetsov et al. TT9^ . GRS 1758-258 also 
does not have yet an identified counterpart, so its mass and or¬ 
bital periods are still unknown. The protoMIRAX observations 
will in principle be capable to attain spectral and flux informa¬ 
tion about GRS 1758-258. 

The GC field to be observed by protoMIRAX will also in¬ 
clude GX l-r4, the best studied accreting pulsar around the GC, 
with a pulse period of the order of 2 minutes. This object is the 
prototype of the small but growing subclass of slow-rotating ac¬ 
creting X-ray pulsars called Symbiotic X-ray Binaries (SyXB), 
in which a neutron star accretes from the wind of a M-type giant 
companion ( Masetti et al.|2006) . GX l-r4 is a persistent source, 
but with strong, irregular flux variations on various timescales 
and extended high/low states. I ts peculiar spin histo ry has been 
the object of intensive study ( [Nagase[[l989[ ), and [Makishima| 
et al.[ ( [1988[ ) have shown a clear transition from spin-up to 


spin-down behaviour at the equilibrium period, characterizing 


a torque re versal episode. More recent studies ( [Gonzalez-Galan 
et al. 2012] ) have used data from BeppoSAX, INTEGRAL, Fermi 


and Swi/Z/BAT to show that the source continues its spin-down 
trend with a constant change in frequency, and the pulse period 
has increased by ~50% over the last three decades. Since the 
source is highly variable, our spectral and timing observations 
can possibly contribute to the flux and period histories for this 
peculiar object. 

Besides providing new data on astrophysical sources, the ex¬ 
periment will be able to make precise measurements of the hard 
X-ray background at balloon altitudes for the Brazilian mid¬ 
latitudes, which is important for investigations concerning the 
South Atlantic Anomaly (SAA) and Space Weather studies. 


3. The protoMIRAX experiment 

The protoMIRAX experiment consists basically in an hard X- 
ray imaging camera mounted in a stratospheric balloon gondola 
capable of carrying out pointing observations of selected target 
regions. Table [^presents an overview of the experiment baseline 
parameters. 

3.1. The X-Ray Camera 

The X-ray camera (XRC) is a coded-aperture wide-fleld hard 
X-ray imager that consists in an aluminum tower-like structure 
that supports an array of X-ray detectors, a collimator, a coded 
mask, passive shields and other electrical shielding and support¬ 
ing components. In this section we describe in details each com¬ 
ponent of the XRC. 


3.1.1. The X-ray detectors 

The X-ray detectors we use in this experiment are band-gap, 
room-temperature semiconductors made of an alloy of 90% 
Cadmium Telluride (CdTe) and 10% Zinc Telluride (ZnTe), 
called CdZnTe or simply CZT. They have high photoelectric ef¬ 
ficiency up to hundreds of keV, due to the high atomic numbers 
of Cd (Z = 48) and Te (Z = 52) and a density of 6 g cm“^. 
The probability of photoelectric absorption per unit pathlength 
is roughly a factor of 4 to 5 times higher than in Ge for typi¬ 
cal gamma-ray energies. The relatively large band gap of 1.5 eV 
precludes significant thermal excitations at room temperature 
and provides good enough energy resolution (AE/E < 10% at 
60keV). 

CZT detectors have been often used in X-ray astronomy due 
to their high efficiency with small thickness (thus reducing back¬ 
ground, which scales with volume) and relative ease of handling 
and mounting, allowing for tiling so that large nearly-contiguous 
detector planes can be built. 

The experiment uses 169 CZT detectors with platinum pla¬ 
nar contacts and dimensions of 10 mm x 10 mm, thickness of 
2 mm each, out of 200 units acquired from eV Products, USA. 
They will be configured in a 13 x 13 array, thus the total area 
of the detector plane is 169 cm^. Due to physical mounting re¬ 
strictions, the edges of adjacent detectors will be separated by 
10 mm. The operational energy range will be from 30 to 200 keV, 
the lower limit determined by residual atmospheric absorption 
at balloon altitudes (~42 km) and the higher limit determined by 
detector photoelectric efficiency given its thickness. 

Each detector is connected directly to a printed circuit 
board (PCB) that carries the front-end analog electronics with 
a preamp, a low-noise amplifier (ENA) and a shaper. The PCBs 
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Table 1. protoMIRAX baseline parameters 


Detector System 


Detector type: CdZnTe (CZT) 
with platinum planar contacts 
Number of detectors: 169 (13 x 13) 
Energy Range: 30 - 200 keV 

Effective area @ 80 keV: 52 cm^ 
(through mask, at 2.7 g cm"^) 

Dimensions: 10 mm X 10 mm 

X 2 mm (thickness) 

Gap between detectors: 10 mm 
Geometrical area: 169 cm^ 

Time resolution: 5 ps 

Collimator 

Blade material: Cu (0.5 mm), 

Pb (0.5 mm), Cu (0.5 mm) 

Blade height: 81 mm 

Cell size: squares of 20-mm sides 

External envelope: 260 mm 

X 260 mm x 81 mm 


Coded Mask 


Material: 1 mm-thick lead 

Basic pattern: 13x13 MURA 

Extended pattern: 2x2 basic 
(minus 1 line and 1 column) 
Position: 650 mm from detector plane 

Open fraction: 0.497 

Element Size: 20 mm x 20 mm 

Total mask dimensions: 500 mm x 

500 mm x 3 mm (total thickness) 
(includes 2mm-thick acrylic substrate) 

Shi 

elding 

Material: lead (external - 1.5 mm), 
copper (internal - 0.5 mm) 

Position: around detectors and 
collimator (sides and bottom) 


Imaging Parameters 


Angular Resolution: 1°43' Total (fully-coded) FOV: 21° x 21° 

Source Location Accuracy: 10'(10 cr) (14.1° x 14.1° FWHM) 

Balloon gondola and flight 


Gondola total mass: ~ 600 kg 
Flight altitude: 42 km (2.7 g cm"^) 
at -23° latitude over Brazil 
Pointing accuracy: ~10 arcminutes 


Gondola dimensions: 1.4 m x 1.4 m 
X 1.8 m (height) 

Single flight duration: < 40 hours 
Star tracker precision: 10 arcseconds 


were designed in our lab and the first 10 were built in house. 
A contractor will replicate them for the whole detector plane. 
Figure shows a photo of one detector and its associated board 
as well as an exploded view of the detector system. The detec¬ 
tors are mounted at 45° with respect to the boards in order to 
minimize the spacing between adjacent detectors. 

The readout circuit of each detector has an applied bias volt¬ 
age of-220 V, a 100 MQ bias resistor and a 150 nF coupling 

capacitor. The charge released by an X-ray interaction in each 
detector will feed an AC-coupled charge-sensitive preamplifier 
that creates an output pulse of typically < 1 mV. The pulses are 
then amplified by ~65 dB by the LNA and formatted by a shaper. 
A Wilkinson-type 8-bit analog-to-digital converter (ADC) con¬ 
verts the pulse amplitudes to time and then measures this time. 
The digitised times are proportional to the energy deposited on 
the detector. This conversion technique is highly linear and pre¬ 
cise; we find it is very suitable for CZT detectors with not very 
high count rates. 

Figure 1^ shows an energy spectrum of a radioactive ^"^^Am 
source for one detector in the lab. The energy resolution 
at the 59.5 keV line is 6.6 keV FWHM (“Full Width at Half 
Maximum”), corresponding to AEjE = 11%. Also shown is 
a peak from a pulse generator that we have used to measure 
the pulse height spread due to electronic noise. The equivalent 
electronic noise resolution is 4.5 keV FWHM. Since the elec¬ 
tronic noise and the Poissonian fluctuations due to the number 
of charge carriers (electron-hole pairs) within the material must 
add in quadrature, the intrinsic, purely statistical resolution is 
4.8 keV FWHM (8%). The red wing on the 60keV peak is due 
to incomplete charge collection within the CZT. The spectrum 
also shows a blend of Cd and Te escape peaks between the ra¬ 
dioactive lines. 


3.1.2. The scientific data acquisition subsystem 


The pulses generated by X-ray interactions in the detectors are 
processed individually by a data acquisition scheme that in¬ 
cludes the front-end analog electronics, described above, and a 
digital electronics for processing and formatting. The detector 
pulses are digitised, tagged with position and time, formatted, 
and then sent to the on-board data-handling computer for stor¬ 
age and transmission to the ground. 

The detector signals of each row of 13 detectors in the detec¬ 
tor plane array feed into a box with the conversion electronics 
(CE) with a Complex Programmable Logic Device (CPLD) and 
microcontrollers. Each one of the 13 CEs acquires data coming 
from the 13 detector PCBs on that particular row and feeds the 
pulses into 13 8-bit ADCs (mentioned above), one for each de¬ 
tector, so the energy of each event is determined in 256 channels. 
The CEs identify the specific detector that is hit on its associated 
row. Each event will be time-tagged unambiguously with 5 yus 
precision by an internal clock whose accuracy is checked each 
second by a GPS receiver. Each CE box receives PPS signals 
from the GPS. 


The signals processed by the CEs will be read by a multi¬ 
plexing electronics (MUX) box at a rate of 1 Hz. The MUX is 
thus responsible for the identification of the detector in which 
the event occurred (the event position), the energy and the rela¬ 
tive time of the event. Eor each acquisition cycle (1 second), the 
MUX will attach to the data package the absolute time informa¬ 
tion coming from the GPS, so each event will have its instant of 
occurrence completely determined within 5 ps. 

The MUX will s end to the on-board data handling computer 
(described in section 3.2.1 ) complete datasets containing all data 
originated in the 13 CEs at a maximum rate of 115,200 bps. In 
this way, the time, position and energy of each event will be 
stored individually so that the position and energy distribution of 
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Fig. 1. Top: one CZT detector (20 mmx20 mmx2 mm) and its as¬ 
sociated electronics board with preamp, LNA and shaper. That 
actual componentes of the LNA are in the opposite side of the 
board to avoid electrical interference. The thick wire is the high- 
voltage supply to the detector (the blue wire is a test lead). 
Bottom: a computer rendering of an exploded view of the de¬ 
tector system. The upper plate is a structural Al frame in which 
a 0.016mm-thick aluminum foil will be stretched to provide light 
and electrical shielding to the detectors. The bottom brass piece 
provides structure and casing for the detectors and PCBs. 



Energy [keV] 

Fig. 2. Energy spectrum of a radioactive Am source for one of 
the CZT detectors. The peak on the right (~ 85 keV) is produced 
by a pulser and has a FWHM of 4.5 keV. The 59.5 keV peak has 
a resolution of 6.6 keV FWHM {^E|E = 11 %). 


events for any desired integration time can be determined. With 
an estimated count rate of ~ 55 counts/s for the entire detector 
plane (see section |^, the dead time will be negligible even for 
very strong sources at levels of several Crabs. 

Figure shows a simplified block diagram of the detector 
data acquisition subsystem. 



Fig. 3. Simplified block diagram of the detector data acquisition 

subsystem. The bias voltage is-200 V, the bias resistor is of 

lOOMQ and the coupling capacitance is of 150 nF. LNA: low- 
noise amplifier; ADC: analog-to-digital converter; MUX: multi¬ 
plexing electronics; OBDH: on-board data handing computer. 


3.1.3. The collimator 

A collimator will be mounted in front of the detector plane in 
order to define the field-of-view (FOV), minimize background 
and provide a uniform illumination fraction of the detectors at 
any incidence angle. In this way, the entire FOV, all the way to 
zero intensity, will be fully coded by the mask, ensuring that 
the experiment does not have any partially coded field-of-view 
(PCFOV). Other coded-mask systems do not use collimators, 
only shielding materials on the sides of the cameras. In these 
systems, sources in the PCFOV may create severe difficulties in 
the reconstruction process by casting incomplete shadowgrams 
of the mask pattern onto the detector plane. In addition, ambigu¬ 
ities in source positions are introduced. We have decided to use 
a collimator to avoid theses issues and reduce background. The 
price to be paid is the reduction of the sensitivity with increas¬ 
ing off-axis incidence angle, according to the collimator angular 
response. 

The collimator cell centers will be separated by 20 mm, in 
such a way that the individual CZT detectors will each be in the 
bottom of a collimator cell. The collimator blades are 81-mm 
high and have lead cores (0.5 mm) and copper plates (0.5 mm) 
on each side to absorb the radiation impinging at high angles 
and also provide a graded shield for lead X-ray fluorescence. 

The collimator is surrounded by a graded shield made of 
lead (1.5 mm) in the outside and copper (0.5 mm) in the inside 
to suppress radiation from outside the field-of-view. The shield¬ 
ing walls have the same height as the collimator blades (81 mm) 
and are separated from the last blades by 2 mm, so that the envi¬ 
ronment that the last detectors (edges and corners) are subjected 
to are very similar to the more central ones. This is important 
to minimize inhomogeneities in the count distribution over the 
detector plane. 


3.1.4. The coded mask 


A coded mask, shown in Figure]^ will be mounted at a distance 
of 650 mm from the detector plane. The mask has the function of 
providing a spatial encoding of the X-ray incoming fluxes (see 
|Dicke|19^ for an introduction of the concept). With the knowl¬ 
edge of the mask pattern and the number of X rays detected by 
each CZT detector for a given integration time, one can recon¬ 
struct an image of the observed field-of-vi ew by suitable decon- 
volution or cross-correlation algorithms (Fenimore & Cannon 
|1978| ). The mask is made out of 1mm-thick lead sheets and the 
basic cells are 2mmx2mm squares. The pattern of closed and 
open elements are a cyclic extension (2x2) of a 13 x 13 Modified 
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Uniformly Red undant Array (MURA) pattern ( [Gottesman & 
|Fenimore|1989 ) with 20mm-side square basic cells. Of the 169 
elements of the pattern, 84 are open and 85 are closed, giving an 
open fraction of 0.497. In order not to have complete repetitions 
of the pattern, which would introduce ambiguities in source po¬ 
sitions in the sky, we remove one row and one column from the 
total 2x2 array of basic patterns, so that the final dimensions of 
the mask are 25 x 25 elements or 500 mm x 500 mm. 



Fig. 4. Photo of the protoMIRAX coded mask. The basic ele¬ 
ments are Imm-thick 20mmx20mm lead squares. The basic 
pattern is a 13x13 MURA repeated four times (2x2) minus 1 
row and 1 column to provide unique shadowgrams for any inci¬ 
dent direction (skybin). The total size is 500mmx500mm. The 
lead pieces are glued to an X-ray transparent, Imm-thick acrylic 
substrate both above and below the mask. 


The MURA patterns belong to a class of aperture patterns 
that allow images to be produced with no intrinsic noise, which 
means that the autocorrelation function of the mask pattern (ac¬ 
tually the cross-correlation between the mask pattern and the de¬ 
coding function, which is almost exactly equal to the mask pat¬ 
tern, with the exception of one single element) is a delta function. 
For a point source, the reconstructed image is “perfect” in the 
sense th at it produces a peak in an other wise completely fiat im- 
age (see [Gottesman & Fenimore| ( |l989| ) and |Braga et al.| ( [200^ 
for details about imaging with MURAs). With this configuration, 
the protoMIRAX X-ray camera will have a total field-of-view of 
20°.8 X 20°.8 totally coded by the mask and a geometrical angu¬ 
lar resolution (one skybin) of 1°43'. Due to the fact that there are 
10-mm gaps between adjacent detectors, the fraction of the FOV 
with maximum sensitivity (total detector area illuminated) will 
be 4°05' X 4°05'. In other words, the collimator blades will only 
start shadowing the detectors for incident directions outside this 
angular range, centered on the camera axis. The source location 
accuracy (see |Caroli et al.|19F7] ) will depend on the incident an¬ 
gle and will be limited by the precision of the pointing system. 
In the best case, it will be around 10' for a lO-cr source in the 
center of the field of view. 


3.1.5. The passive shield 

The detector array and associated electronics, including the de¬ 
tector boards, power supplies and batteries, will be surrounded at 
the sides and bottom by a gradual passive shield of lead (1.5 mm) 
and copper (0.5 mm). This shield will be responsible for absorb¬ 
ing most of the radiation coming from directions other than those 
defined by the collimator. Figure shows computer diagrams 
with the main parts of the XRC, including the aluminum support 
structures. 



Fig. 5. Computer rendering of the protoMIRAX X-ray Camera. 
Left: Exploded diagram of the XRC, showing, from top to bot¬ 
tom, the coded mask, collimator, aluminum foil support plate, 
detector system, detector electronics and power supply boxes, 
and battery pack. Right: The assembled XRC with passive 
shielding and structural parts, including the elevation axis sup¬ 
port. Also shown are one of protoMIRAX’s two star cameras 
(top right) and the elevation driving mechanism (bottom left) 
based on a ball screw. 


3.2. The Balloon Gondola 

The protoMIRAX balloon gondola will house the X-ray cam¬ 
era and the various subsystems of the space segment, in¬ 
cluding the On-Board Data Handling Subsystem (OBDH), the 
Attitude Control and Pointing Subsystem (ACS), the Telemetry 
and Command Subsystem (TM&TC), and the Power Supply 
Subsystem (PSS). Also, there will be a main GPS receiver that 
will provide universal time to all subsystems and two star cam¬ 
eras that will be part of the control loop. A computer rendering 
of the gondola is shown in Figure 

The gondola will use a power supply system based on pack¬ 
ages of Li-ion batteries. A telemetry system operating in the 
L-band (1.5 GHz) will provide data and housekeeping links to 
the ground. A command system will provide pointing opera¬ 
tions during the balloon flight. The science data will be stored 
on-board on computer flash memories and will also be entirely 
transmitted to the ground. 

3.2.1. The On-Board Data Handling Subsystem 

The On-Board Data Handling Subsystem (OBDH) is responsi¬ 
ble for acquiring, formatting and transmitting all data that come 
from the various subsystems of protoMIRAX’s space segment 
to the ground station (GS). The OBDH is also responsible for 
receiving and retransmitting, when necessary, the various com¬ 
mands sent by the GS to the gondola. For each X-ray photon 
detected by the XRC, a 6-Byte packet is created encasing the 
time stamp, x-y position (detector that was hit) and energy (pulse 
height) of the event that were formatted and sent by the MUX. 
These event packets are then sent every 1 second to the OBDH 
through a 115.2 kbps RS-422 unidirectional serial communica¬ 
tion line. In addition, another 115.2 kbps RS-422 serial interface 
allows OBDH to send specific commands to the XRC and re¬ 
ceive XRC’s housekeeping data. 

We devised a cost-effective and robust solution for the ex¬ 
periment timing in which a single Global Positioning System 
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Fig. 6. Computer rendering of the protoMIRAX balloon gon¬ 
dola, showing the X-ray camera in the front face. The coloured 
boxes at the bottom of the frame are battery packs and electron¬ 
ics racks. The flat box in the back is the TM&TC subsystem. The 
reaction wheel can be seen below the XRC. One can also see the 
two star cameras: one is fixed at the XRC and the other (behind 
the coded mask) is fixed in the gondola frame. At the top it is 
shown the decoupling mechanism that connects to the balloon 
cables. The four brown boxes in the bottom are crash pads. The 
gondola frame dimensions are 1.4 m x 1.4 m x 2 m. 


(GPSDXA) unit sends PPS (1 pulse per second) and data to the 
various subsystems of the experiment: OBDH, XRC, ACS, and 
the Autonomous Star Trackers (AST). The OBDH attaches the 
UT data from the GPSDXA to the scientific data files at every 
second, so absolute timing information at the resolution of 5 yus 
can be retrieved for every event. 

The OBDH also communicates with the GS and the 
ACS. A 500kbps synchronous communication channel connects 
the OBDH with the Flight Control and Telecommunications 
Subsystem (FCTS). This channel is exclusively dedicated to 
transmit scientific data to ground. FCTS also communicates with 
the GS via a double transmitter/receiver system operating in L- 
band. In order to transmit all housekeeping data generated by 
the several subsystems of the space segment, OBDH uses a se¬ 
rial RS-232 channel operating at 115.2 kbps. Ground station 
commands are received by OBDH via FCTS by means of a 9.6 
kbps RS-232 serial channel. In addition, all ACS data are sent to 
OBDH prior to be sent to the GS via another 115.2 kbps RS-422 
serial interface standard. 

The main hardware unit of OBDH is the Payload Data 
Handling Computer (PDCpM), which is a PC/104 ultra low 
power AMD Geode LX computer. It operates on 333 MHz, has a 
128-MByte SDRAM, and has all the necessary interfaces previ¬ 
ously mentioned in addition to analog-digital/digital-analog con¬ 
verters. The PDCpM went through thermal cycling and thermal 
vacuum testing with temperatures ranging from -40° C to -r85° 
C and pressure as low as 1.3 mbar. 


The software embedded into the computer (SWPDCpM) has 
been developed in C language over the real-time operating sys¬ 
tem RTEMS. The architecture of the SWPDCpM is composed 
of 3 layers where the bottom layer is the basic software, the in¬ 
termediate layer is the flight software library in which there are 
basic services and reusable components, and the top layer is the 
application software, the main part of SWPDCpM. 

Figure shows a simplified physical architecture of the pro¬ 
toMIRAX experiment with its main subsystems. 
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Fig. 7. Simplified physical architecture of the protoMIRAX 
balloon experiment. XRC: X-ray camera; ACS: Attitude 
Control System; TEMPDXA: Temperature monitoring equip¬ 
ment; GPSDXA: GPS unit; OBDH: On-Board Data Handling 
Subsystem; PDCpM: Payload Data Handling Computer; DC- 
DC Conv: DC-DC Converter; PRESN: Pressure Sensor; 
PSS: Power Supply Subsystem; FCTS: Flight Control and 
Telecommunications Subsystem. 


3.2.2. The Attitude Control and Pointing Subsystem 

The X-ray camera will be placed in an alt-azimuth mount in the 
balloon gondola, allowing pointing to different targets in the sky 
and tracking them continuously down to elevation angles of 30°. 
A ball screw mechanism, with a step motor, will allow smooth 
and secure motion in elevation, and the azimuth motion will be 
provided by a system that rotates the entire gondola with respect 
to the geomagnetic field and the positions of the sun and the stars 
in the sky. The attitude control system will employ several sen¬ 
sors and actuators. The sensors include an electronic compass, 
an accelerometer, two star sensors, a sun tracker and both az¬ 
imuth and elevation angular encoders. The actuators include a 
reaction wheel for angular momentum exchange and dump, and 
motors that will drive the motions in both elevation and azimuth. 
The gondola also has a mechanism at its top for decoupling its 
rotation with respect to the balloon itself and a motor for desat- 
urating the reaction wheel. The overall pointing accuracy of the 
control system is < 0.5°, depending on the elevation angle; high 
elevations are in general more difficult to control since small 
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changes in the pointing direction near the zenith leads to large 
azimuthal and roll angle displacements. Also, the accuracy will 
depend on whether the flight is at daytime or nighttime. At night¬ 
time, the star cameras will provide a more accurate control. 

The two autonomous star sensors of protoMIRAX are be¬ 
ing developed at INPE. They are wide-field (25°x25°) optical 
(400-800 nm) cameras with a 35 mm focal length. The imaging 
sensor is a CMOS APS (active pixel sensor) with 1024x1024 
pixels and quantum efficiency up to 25%. An especially de¬ 
signed algorithm for star pattern recognition provides attitude 
knowledge precision better than lOarcseconds (1 cr) around the 
cross-foresight axes (yaw and pitch) and better than 60 arc sec¬ 
onds (1 cr) around the foresight axis (roll angle). One of the star 
sensors is fixed in the XRC mechanical frame so that it enables 
precise pointing of the XRC telescope. Its optical axis is tilted 
24° with respect to the X-ray axis, toward lower elevations, so 
as to avoid having the balloon in the field of view for high ele¬ 
vation observations. Even with this misalignment, there will be 
occasions when the primary star tracker will be obstructed by the 
balloon. In these cases the XRC attitude will be determined by 
the secondary star tracker mounted in the gondola, albeit with a 
somewhat reduced accuracy due to uncertainties from the angu¬ 
lar encoders used in the motors that control the XRC mechanical 
frame elevation. The star trackers are also capable of acquiring 
images, which will be used for evaluation of the star tracker per¬ 
formance during the flight. 


3.3. The Ground Subsystem 

The operations of the protoMIRAX telescope during the balloon 
flights will make use of the ground facilities already available 
at INPE. The reuse strategy focused on two subsystems: (i) the 
balloon flight control (OPS) and (ii) the experiment operation it¬ 
self, which will be supported by the SATellite Control System 
(SATCS). SACTS is a software-based architecture developed 
at INPE for satellite commanding and monitoring. Considering 
that protoMIRAX is a pathfinder for the MIRAX satellite mis¬ 
sion, a ground infrastructure compatible with INPE’s satellite 
operation approach would be useful and highly recommended 
to control and monitor the experiment during the balloon cam¬ 
paigns. 

The OPS required investments on information technology. 
The available balloon ground station needed little hardware 
maintenance, but significant efforts on software updating were 
necessary. The original standalone computer system dedicated to 
run the flight control software was modernized both in hardware 
and software by the new system named OPS/ES. In addition, 
a new server computer, properly configured for Ethernet con¬ 
nections, has extended the existing ground station facilities with 
a network switch, serial converters and a new software named 
OPS/Server in order to support the available uplink and down¬ 
link channels being mapped to TCP/IP gateways. Those com¬ 
munication improvements were necessary to support the inter¬ 
operability between the balloon ground station, SATCS and pro¬ 
toMIRAX data center. 

Concerning SATCS, some effort in software customization 
was necessary because SATCS architecture complies with par¬ 
ticular operational requirements on different missions by us¬ 
ing several customized object-oriented software elements and 
frameworks. The diagram in Eigurej^ shows the ground solution 
designed for the Ground Segment of the protoMIRAX experi¬ 
ment. 
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Fig. 8. Physical architecture of the protoMIRAX ground subsys¬ 
tem. 

4. Effective Area and Sensitivity 

The effective area of a balloon experiment for the observation of 
cosmic sources is a function of energy (E), atmospheric depth 
and zenith angle: 

( 1 ) 

where Ageo is the geometrical detection area, 6=1- exp(-/jD 1) 
is the detector efficiency, fioiE) is the detector’s attenuation co¬ 
efficient in cm“\ I is the detector thickness in cm, {iiIp)a{E) is 
the attenuation coefficient of the atmosphere in cm^/g, x is the 
atmospheric depth in g cm“^ and z is the zenith angle. 

Eor the expected observations of protoMIRAX at 2.7 g cm“^, 
considering a point source at the zenith, the effective area is 
shown in Eigurej^ with a maximum value of 52 cm^ at ~80keV. 
It is important to note that the geometrical area for a point source 
at infinity is the total detector area of 169 cm^ multip lied by the 
open fraction of the mask, which is 0.497 (see sectior (3X4| ). 



Energy (keV) 


Fig. 9. Effective area of protoMIRAX. 


The sensitivity of the experiment can be calculated by esti¬ 
mating the background at the altitudes expected for the balloon 
flight. This was produced by performing Monte Carlo simula¬ 
tions using the software package GEANT4 ( [Agostinelli et al 


|2003| ), a well-known suite of routines that performs detailed cal¬ 
culations of gamma-ray and particle interactions on all materials. 
Given the experiment mass model and the environmental par¬ 
ticle and photon fields, the code provides the detector spectral 
response, both prompt and delayed. 
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Based on detailed GE ANT4 calculations we have carried out 
( Penacchioni et al.| ( |2015| ); Castro et al., in preparation), the esti¬ 
mated background count rate for the entire detector plane, in the 
30 to 200 keV energy range, is of ~ 55 counts/s. In these sim¬ 
ulations, we have considered both the anisotropic atmospheric 
X and y-ray spectra and the cosmic diffuse contribution com¬ 
ing from the telescope aperture. We have also considered pri¬ 
mary and secondary protons, electrons and neutrons. As shown 
by |Penacchioni et aL] ( |2015| ), the photon contribution is dominant 
up to 90 keV, where the contribution due to protons becomes al¬ 
most equally important. Neutrons are important below 35 keV 
and the electron contribution is negligible. 

In Figure 10 we show a map of the simulated background 
counts over the detector plane for an integration of 4 hours at 
2.7 g cm“^ for the entire energy range 30-200 keV. The in¬ 
creased count rates at the edges and corners of the detector plane 
are most likely due to the closeness of these detectors with the 
shielding materials and support structures, which produce sec¬ 
ondary radiation by scattering. Another factor could be due to 
the extra shielding produced by the collimator blades to the de¬ 
tector in the central areas of the detector plane. 


Background 



Fig. 10. Background count distribution over the detection plane 
in the detector energy range 30 - 200 keV, for an observation 
time of 4 h at balloon altitudes. 


Taking the expected background levels at balloon altitudes 
as a function of energy, one can calculate the sensitivity of the 
experiment in terms of the minimum detectable flux at a particu- 
lar stat istical significance. As shown by |Gottesman & Fenimore| 
( 1989| ), the signal-to-noise-ratio (SNR) of a coded mask instru¬ 
ment that uses URAs or MURAs, for a single point source, is 
given by 


SNR= — -^- . (2) 

where Ns is the net source counts and Nb is the background 
number of counts for a given integration time in a particular en¬ 
ergy range. Solving for Ns , one gets 

Ns = SNR {SNR + ^JSNR^+NB). (3) 

If we are interested in a minimum detectable flux, the SNR will 
be a small number (< 10), and it is reasonable to assume that 
SNR^ Nb for sufficiently long integration times. Therefore, 
Ns ~ SNR ^[N^. 


Now, if S is the source flux in photons cm“^s“^keV“\ then 
Ns = S Aeff T AE, where Aeff is the effective area in cm^, T is 
the integration time in seconds and AE is the energy interval un¬ 
der consideration (in keV). Similarly, Nb = BA^qq T AE, where 
B is the background level in the more suitable units of counts 
cm“^s“^ and Ageo is the geometrical area of the whole detector 
plane. 

The minimum source flux that will be detectable at a level of 
No- = S NR will then be 




Aeff AE 


BA 


geo 


photons cm^s ^keV ^ 


(4) 


According to this expression, the 3-cr sensitivity (minimum de¬ 
tectable flux) for the whole energy band is found to be approx¬ 


imately 1.9 X 10“^ counts cm“"^ s“^ from 30 to 200 keV, con¬ 
sidering 8 h of integration at 2.7 g cm“^ and a zenith angle of 
30°. This allows protoMIRAX to make detailed observations of 
the Crab (even at low elevation angles) and also detect a few X- 
ray binaries in states of high hard X-ray emission. In Figure [TT] 
we show Si 3-cr sensitivity curve for protoMIRAX for an integra¬ 
tion time of 8 hours at 2.7 g cm“^, considering a point source 
at a zenith angle of 30°. The increased minimum flux values be¬ 
tween 70 and 90keV are due to somewhat strong lead fluores¬ 
cent lines (reported by |Penacchioni et aL] ( |2015| )) in the estimated 
background spectrum. We are in the process of studying differ¬ 
ent shielding configurations that will very likely minimize this 
problem. 


protoMIRAX on-axis sensitivity 



Fig. 11. On-axis sensitivity curve for protoMIRAX. The hori¬ 
zontal bars are the minimum detectable fluxes at a level of 3cr 
at an atmospheric depth of 2.7 g cm“^ and a zenith angle of 
30°. The integration time is 8 hours. Also shown are the spec¬ 
tra of 4 X-ray sources that will be observed in the first balloon 
flight. The source spectra were taken from |Sizun et ak] (|2004|) 
(Crab), Grebenev et al. ( |I995| ) ( IE 1740.7-2942), |Sunyaev et al. 


( |1991b| ) (GRS 1758 - 258) and |Dieters et"^ ( |1991| ) (GX 1 -r 4). 


5. Image Simulation 

In order to predict the performance of the protoMIRAX imaging 
system, we have simulated an observation of the Crab nebula 
during a meridian transit at our flight latitudes (23° S). To calcu¬ 
late the number of counts from the source, we have considered 
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the zenith angle variation with time before and after the tran¬ 
sit, which occurs at an elevation of 45°. Based on the Crab hard 
X-ray spectrum taken from Sizun etak] ( |2004| ), we built with 
GEANT4 a simulated shadowgram of a Crab observation in the 
center of the FOV, which is shown in Figure 



Fig. 12. A count distribution (shadowgram) on the detector plane 
of a simulated observation of the Crab, including the expected 
background, for an integration time of 4 h around the meridian 
passage at a latitude of 23° S. The energy range is 30 to 200 keV. 


The image covers the entire FOV of 21° x 21°, and each 
sky bin corresponds to 1°43' of the sky. The overall signal-to- 


noise ratio (SNR) of the image is 109 cr (see Penacchioni et al. 
|2015[ for details). This value corresponds to 91% of the theoret¬ 
ical value given by expression (2) above. The loss in detection 
significance is mainly due to the non-uniformity of the back¬ 
ground across the detector plane. In Figure \T3\ we show the re¬ 
constructed Crab image. It is interesting to note that the Crab 
will be in principle observable by protoMIRAX at a 5-(T level 
for an integration as short as 63 s, which makes it a very good 
calibration source for the experiment. 



Fig. 13. Simulated image of the Crab Nebula region for an in¬ 
tegration time of 4 h around the meridian transit at balloon alti¬ 
tudes (2.7 g cm“^). The energy range is 30-200 keV. The signal- 
to-noise ratio of the detection is 109. 


6. Conclusion 

The protoMIRAX experiment is a balloon-borne hard X-ray 
telescope being developed as a scientific and technical pathfinder 


for the MIRAX satellite mission. In the current configuration, 
under revision due to restrictions imposed by the Brazilian space 
program, MIRAX will consist in a set of wide-field coded-mask 
hard X-ray cameras (5-200 keV) with an angular resolution of 
~ 5-6 arcmin that will operate in scanning mode in a near- 
equatorial circular low-Earth orbit. MIRAX main goals include 
the characterization, with unprecedented depth and time cover¬ 
age, of a large sample of variable and transient phenomena on ac¬ 
creting neutron stars and black holes, as well as Active Galactic 
Nuclei and GRBs. 

Apart from serving as a testbed for several MIRAX subsys¬ 
tems, protoMIRAX will be capable of doing interesting science 
of its own. The hard X-ray camera, similar to the hard X-ray 
imager being developed for MIRAX, albeit with a much lower 
angular resolution and smaller detector area, has a resolution of 
1°43' with a fully-coded FOV of 21° x21°. Detailed background 
calculations have been performed with the GEANT4 package, 
and the expected background at balloon altitudes at a latitude of 
-23° over the SAA region in Brazil provides a 3-cr sensitivity 
for the 30-200 keV range of ~ 1.9 x 10“^ photons cm“2 s“^ 
for an integration of 8 hours. In its first balloon flight, protoMI¬ 
RAX will observe the Crab Nebula for calibration and imaging 
demonstration and also the GC region for scientific purposes. 
The experiment is in a very advanced state of development and 
we hope to carry out the first ballon flight in late 2015. 

Acknowledgements. We thank FINEP, CNPq and FAPESP, Brazil, for financial 
support. We also thank Fernando G. Blanco, Sergio Admirabile, Luiz A. Reitano, 
Wendell P. da Silva, Leonardo Pinheiro, Fernando Orsatti, Paulino Scherer and 
the engineers from the company COMPSIS in Sao Jose dos Campos for in¬ 
valuable technical support. A.V.P. acknowledges the support by the interna¬ 
tional Cooperation Program CAPES-ICRANET financed by CAPES - Brazilian 
Federal Agency for Support and Evaluation of Graduate Education within the 
Ministry of Education of Brazil. M.C.A. and J.R.S. also acknowledge CAPES 
for Ph.D. and post-doctoral fellowships, respectively. 


References 

Agostinelli, S., Allison, J., Amako, K., et al. 2003, Nuclear Instruments and 
Methods in Physics Research A, 506, 250 
Braga, J., D’Amico, K, Villela, T, et al. 2002, Review of Scientific Instruments, 
73,3619 

Braga, J., Rothschild, R., Heise, J., et al. 2004, Advances in Space Research, 34, 
2657 

Biihler, R. & Blandford, R. 2014, Reports on Progress in Physics, 77, 066901 
Caroli, E., Stephen, J. B., Di Cocco, G., Natalucci, L., & Spizzichino, A. 1987, 
Space Sci. Rev., 45, 349 

Castro, M., D’Amico, K, Braga, J., et al. 2014, A&A, 569, A82 
Dicke, R. H. 1968, ApJ, 153, LlOl 

Dieters, S. W. B., Greenhill, J. G., Sharma, D. P, et al. 1991, Advances in Space 
Research, 11,35 

Fenimore, E. E. & Cannon, T. M. 1978, Appl. Opt., 17, 337 
Gehrels, N. & Swift Team. 2004, New A Rev., 48, 431 
Gonzalez-Galan, A., Kuulkers, E., Kretschmar, P, et al. 2012, A&A, 537, A66 
Gottesman, S. R. & Fenimore, E. E. 1989, Applied Optics, 28, 4344 
Grebenev, S. A., Pavlinsky, M. N., & Sunyaev, R. A. 1995, Advances in Space 
Research, 15, 115 

Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013, ApJ, 770, 103 
Jager, R., Mels, W. A., Brinkman, A. C., et al. 1997, A&AS, 125, 557 
Kuznetsov, S., Gilfanov, M., Churazov, E., et al. 1996, in Roentgenstrahlung 
from the Universe, ed. H. U. Zimmermann, J. Triimper, & H. Yorke, 157-158 
Levine, A. M., Bradt, H., Cui, W., et al. 1996, ApJ, 469, L33 
Liang, E. P. & Nolan, P. L. 1984, Space Sci. Rev., 38, 353 
Madsen, K. K., Reynolds, S., Harrison, F., et al. 2015, ApJ, 801, 66 
Makishima, K., Ohashi, T, Sakao, T, et al. 1988, Nature, 333, 746 
Masetti, N., Orlandini, M., Palazzi, E., Amati, L., & Frontera, F. 2006, A&A, 
453, 295 

Mirabel, I. F., Rodriguez, L. F, Cordier, B., Paul, J., & Lebrun, E. 1992, Nature, 
358,215 

Nagase, F. 1989, PASJ, 41, 1 

Paul, J., Ballet, J., Cantin, M., et al. 1991, Advances in Space Research, 11, 289 


9 








































J. Braga et aL: The protoMIRAX balloon experiment 


Penacchioni, A. V., Braga, J., Castro, M. A., & D'Amico, R 2015, Journal of 
High Energy Astrophysics, in press 

Ricker, G. R., Atteia, J.-L., Crew, G. B., et al. 2003, in American Institute 
of Physics Conference Series, Vol. 662, Gamma-Ray Burst and Afterglow 
Astronomy 2001: A Workshop Celebrating the First Year of the HETE 
Mission, ed. G. R. Ricker & R. K. Vanderspek, 3-16 
Sizun, P, Shrader, C. R., Attie, D., et al. 2004, in ESA Special Publication, 
Vol. 552, 5th INTEGRAL Workshop on the INTEGRAL Universe, ed. 
V. Schoenfelder, G. Lichti, & C. Winkler, 815 
Sunyaev, R., Churazov, E., Gilfanov, M., et al. 1991a, ApJ, 383, L49 
Sunyaev, R., Gilfanov, M., Churazov, E., et al. 1991b, Soviet Astronomy Letters, 
17, 50 

Sunyaev, R. A., Babichenko, S. L, Goganov, D. A., Tabaldyev, S. R., & 
lamburenko, N. S. 1990, Advances in Space Research, 10, 233 
Takahashi, T, Mitsuda, K., Kelley, R., et al. 2014, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 9144, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 25 
Tavani, M., Bulgarelli, A., Vittorini, V, et al. 2011, Science, 331, 736 
Virgilli, E., Frontera, F., Valsan, V, et al. 2013, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 8861, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 7 
Wilson-Hodge, C. A., Cherry, M. L., Case, G. L., et al. 2011, ApJ, 727, L40 
Winkler, C. 1995, Experimental Astronomy, 6, 71 


10 



